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We study three basic diffusion-controlled reaction processes—annihilation, coalescence, and ag¬ 
gregation. We examine the evolution starting with the most natural inhomogeneous initial configu¬ 
ration where a half-line is uniformly filled by particles, while the complementary half-line is empty. 
We show that the total number of particles that infiltrate the initially empty half-line is finite and 
has a stationary distribution. We determine the evolution of the average density from which we 
derive the average total number N of particles in the initially empty half-line; e.g., for annihilation 
{N) = ^ + For the coalescence process, we devise a procedure that in principle allows one to 
compute P{N), the probability to find exactly N particles in the initially empty half-line; we com¬ 
plete the calculations in the first non-trivial case {N =1). As a by-product we derive the distance 
distribution between the two leading particles. 

PACS numbers: 05.40.-a. 82.20.-w. 66.10.C—. 05.70.Ln 


I. INTRODUCTION 

Reaction-diffusion systems are ubiquitous in biology, 
chemistry and physics. An important class of such sys¬ 
tems, diffusion-controlled processes in which a reaction 
happens whenever two reactants “meet”, has been stud¬ 
ied for almost a century (see e.g. [THI3]). Mathemat¬ 
ically, these processes are strongly interacting infinite 
particle systems, so it is not surprising that many ba¬ 
sic questions remain unanswered. Here we analyze a few 
basic diffusion-controlled processes focusing on the inter¬ 
play between reaction and spatial inhomogeneity of the 
initial setting. We will see that for an extreme inhomo¬ 
geneous initial configuration with an empty half-line Q , 
the statistics of the occupation number of this half-line 
is asymptotically stationary and highly non-trivial. 

We study three diffusion-controlled reaction systems. 
One is a single-species annihilation which is represented 
by the reaction scheme 

A + A^lb ( 1 ) 

Identical particles (denoted by A) undergo diffusion and 
whenever two particles touch each other, they disappear. 
Strict annihilation is rare. One example is the annihila¬ 
tion of domain walls in an Ising spin chain subjected to a 
zero-temperature spin-flip dynamics m- In most situa¬ 
tions, however, annihilation merely implies that the reac¬ 
tion product does not further affect the reaction process. 
For instance, a collision between two atoms may lead to 
the formation of a diatomic molecule; if such molecules 
are stable (that is, they do not break back into atoms) 
and if they do not influence diffusing atoms, we can use 
the reaction scheme 0. 

Another simple diffusion-controlled process, coales¬ 
cence, can be represented by the reaction scheme 

( 2 ) 

Strict coalescence is also rare. It occurs for example in 
a Potts chain with infinitely many states evolving ac- 
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FIG. 1: Illustration of the irreversible reaction process. The 
initial condition with fixed concentration of particles left of 
the origin and no particles to the right of the origin is shown 
on the top. The next two layers illustrate the state of the 
system at later times (the latest time at the bottom) with a 
dwindling concentration of particles left of the origin and a 
fluctuating finite total number of particles to the right of the 
origin. 


cording to the non-conserved zero-temperature dynam¬ 
ics [13]. In most situations, the coalescence process is a 
coarse-grained description of the aggregation process 

Ai -h Aj —>■ Ai^j (3) 

where Am represents clusters of mass m, i.e., clusters 
composed of m monomers. The symbolic representation 
([^ indicates that we disregard the mass of aggregates 
and hence it faithfully describes the gross features of the 
aggregation process when the diffusion constant is mass- 
independent and when the size is also mass-independent. 

Little is known about evolution of irreversible reac¬ 
tions starting from spatially inhomogeneous configura¬ 
tions. Here we consider the simplest situation when at 
time t = 0 the half space a; < 0 is uniformly filled by 
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particles, while the half space a; > 0 is empty (see Fig. 


c{x, t = 0) 


p X <Q 

0 a: > 0 


( 4 ) 


This step-function initial condition has been used e.g. to 
study front propagation in reaction-diffusion processes 
[H IS] and shock waves in driven exclusion processes 
[16j . it is also popular in experimental [niiiH] and the¬ 
oretical (see e.g. studies of cold quantum gases, 

and in analyses of quantum spin chains f |23H26j and ref¬ 
erences therein). 

In three dimensions, one can use the rate equation ap¬ 
proach to analyze the evolution. For annihilation and 
coalescence processes, the density c{x,t) satisfies the 
reaction-diffusion equation, dtc = Dd^xC — Kc^ , where D 
is the diffusion coefficient and K the reaction rate. The 
long time behavior is simple, viz. in the scaling limit 

X 

\x\ —>■ oo, t —>■ oo, X = = finite, (5) 

y2Dt 

the density becomes c{x,t) = {Kt)~^C{X). The scaled 
density C{X) satisfies an ordinary differential equation 
which can be solved numerically. 

In this paper we focus on the one-dimensional case. 
This is the most challenging setting and the emergent re¬ 
sults are richer than in higher dimensions. The challenge 
stems from the fact that in one dimension one cannot 
employ the rate equation approach, i.e., one cannot rely 
on reaction-diffusion equations. The behaviors are still 
well understood when the initial conditions are homoge¬ 
neous, e.g., c ^ see e.g. [i El [H Eg. For 

the inhomogeneous initial condition (Fig. the density 
eventually decays in the same manner everywhere. How¬ 
ever, particles infiltrate only up to distance ~ \/~Di into 
the vacant half-line, so the total number of particles in 
the initially vacant half-line is expected to be finite and 
governed by a stationary distribution. In this paper, we 
establish this behavior analytically. 

For the inhomogeneous setting Q, the scaling form 

is anticipated. We shall obtain the scaled density profile 
C{X) analytically. An immediate consequence of is 
that the average total number of particles (N) in the 
initially empty half-line is finite. More precisely, for the 
coalescence process 


(N), = 5 ^ = 0.534154943... (7) 

For the annihilation process, the value is twice smaller. 
Thus the total number of particles N in the initially 
empty half-line is a fluctuating quantity with a finite av¬ 
erage. (The total number of particles in the initially filled 
half-line is of course infinite, although the density van¬ 
ishes algebraically with time ^ Furthermore, 


we also find that the total number of particles N is gov¬ 
erned by a stationary probability distribution P{N). 

The determination Pc{N) for the coalescence process is 
technically challenging. In addition to Pc(0) = | which 
can be established on the basis of symmetry alone, we 
are able to compute the probability 


Pc{l) 
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0.466095976.. 


arctan 



— 2 arctan 



( 8 ) 


to have exactly one particle in the initially empty half¬ 
line. With much more effort it could be possible to de¬ 
termine F’c(2), but to find the entire distribution Pc{N) 
requires more advanced methods. 

For the annihilation process we haven’t been able even 
to determine Pa{0). Indeed, the coalescence process is 
more tractable than the annihilation process. This sub¬ 
tle difference between these two very similar diffusion- 
controlled reaction processes was already encountered in 
the homogeneous setting. For instance, for the coales¬ 
cence process the distance distribution between adjacent 
particles is well-known (see mm). while for the an¬ 
nihilation process the distance distribution has not been 
determined analytically [571 - 12^ . Another long-standing 
unresolved problem is to determine the long-time behav¬ 
ior of an impurity (a particle with diffusion coefficient 
generally different from the diffusion coefficient of the 
host particles) in the annihilation process [301 El] ■ 

The annihilation process with the step-function initial 
condition was investigated in Ref. [3g. That study 
was focused on the survival probabilities 5'n(t), where 
n = 1 , 2 ,... labels particles according to their initial lo¬ 
cations: a;i(0) > X 2 ( 0 ) > 2 : 3 ( 0 ) etc. These survival prob¬ 
abilities exhibit an intriguing temporal behavior, the de¬ 
cay laws depend only on the parity of the original label: 
Sn ^ for even n and as 5'„ ~ t~^ for odd n; the 
decay exponents a and /3 are still unknown [32] . 

The rest of this article is organized as follows. In 
Sec. [n] we consider diffusion-controlled annihilation, coa¬ 
lescence, and aggregation processes in one dimension and 
compute the average densities everywhere. Statistics of 
the total number of particles in the initially empty half¬ 
line is studied in Sec. IIIIl In Sec. IIVI we summarize our 
results and discuss several open problems. 


II. DENSITY PROFILE 


In one dimension, we compute the density profile us¬ 
ing the empty interval method which allows an efficient 
analytical treatment of the one-dimensional diffusion- 
controlled coalescence process (Sec. |II A ). Modifications 
of the empty interval method which are suitable for an¬ 
nihilation and aggregation are presented in Sec. |IIB|and 

sec.rrrci 

We consider continuous versions of the aggregation, 
annihilation and coalescence processes in which size-less 
particles undergo Brownian motions on a line. Lattice 
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versions in which each particle occupies a site on the 
one-dimensional lattice and hops to nearest-neighboring 
sites are also tractable. The basic asymptotic behaviors 
are identical in both versions. 


A. Coalescence Process 


In the one-dimensional diffusion-controlled coalescence 
process, point particles undergo identical independent 
Brownian motions and merge instantaneously whenever 
they meet. An elegant exact treatment of this strongly 
interacting infinite particle system is possible through the 

empty interval technique [101 HU [131 HZl I33H3S] • 

To explain the empty interval technique, we begin with 
spatially-homogeneous coalescence process. Let E{x, y; t) 
be the probability that the interval [x, y] is empty at time 
t. In the homogeneous case these probabilities depend 
only on the length of the interval: E{x,y;t) = E{i,t) 
where £ = y — x. The probabilities E{£, t) satisfy the 
diffusion equation 


dE{e,t) d^E{i,t) 

dt d£^ 

and the boundary condition 


(9) 


E{£ = 0,t) = l (10) 

The initial condition is E{£,t = 0) = if the initial 
distribution of particles is random and the initial density 
is p. Note also the general formula for the concentration 


c{t) 


dE{£, t) 
di 




( 11 ) 


The physically relevant region is ^ > 0, but we can 
extend 0 to ^ < 0. It proves convenient to make a shift, 
F(l, t) = E{1, t) — 1, so that the boundary condition (10) 
becomes F{£ = 0,t) = 0. The governing equation for 
F(£,t) remains the diffusion equation 


OF 


2D 


52 F 
9^ 


( 12 ) 


We have F{£, t = 0) = — 1 when £ > 0 and we seek 

the initial condition for £< 0 in such a way that the 
boundary condition F{£ = 0, t) = 0 is manifestly obeyed. 
The proper choice is 


The solution to (12) and (14) is F = —FTf{£/^/8Di), and 
therefore 


E{£, t) = Erfc 


VWDtJ 


(15) 


By combining (11) and (15) we recover the long-time 
asymptotic behavior of the density 


z(t) = 


1 


y/2'KDt 


(16) 


We now turn to the spatially-inhomogeneous initial 
condition The mathematical formulation is a natural 
generalization of Eqs. 0 -"([TT|). The governing equations 


are 


The boundary condition becomes 
E{x, x]t) = 1 


(17) 


(18) 


The density is found from 


c{x, t) 


dE{x,y;t) 

dy 


y^x 


(19) 


We have to solve 0 subject to the initial condition 


{ 1 0 < X <y 

exp [per] X <Q <y (20) 

exp[p(a: - y)] x < y < 0 


Writing E{x,y;t) = 1-1- F{x,y;t) and focusing on the 
long-time limit we can again use a simplified initial con¬ 
dition: 


Fo{x,y) 


{ 0 0 < min(x, y) 

— 1 X < 0, X < y 
1 y <0, y <x 


Changing variables 

^ = y + x, r] = y-x 


( 21 ) 


( 22 ) 


we end up with the two-dimensional diffusion equation 


F{£,t = 0) = 


exp[—p£] — 1 £ > 0 

1 — exp[p£] £ <0 


(13) 


One can write an exact solution of the diffusion equation 
0 subject to the initial condition ( |I^ . We do not 
display this solution since our main goal is to establish the 
long-time behavior, and for that purpose, we can replace 
(13) by the simplified initial condition 


F{£,t = 0) = 



|-F(C,y;t) = 2F 


52 92 

-b 


F(C,y;t) (23) 


ar \de 5772 

subject to the initial condition [Fig. 

[0 o<^, 

= 1“1 

I 1 p < 0, ^<-r] 


(24) 


Expressing the boundary condition (18) in terms of the 


auxiliary function F(^, rj; t) gives F{^,r] = 0; t) = 0. This 





















4 



FIG. 2: The initial condition fb(C,??) = = 0) in 

the (^, 77 ) plane. The vanishing of F{^, rj; t) on the ^ axis, 
_F|^=o = 0, is manifest. 


boundary condition is indeed manifestly obeyed with the 
choice of initial condition (241, see Fig. 

The solution to equation (23)-([24| is straightforward. 
One gets 


E{x,y;t) = 1 - ^Erfc Erf 


ru+Y 


du e 


dv e 


lu+X 


(25) 


where X = x/^/2Dt and Y = y/VXDt. By substituting 
(251 into (191 we confirm the scaling form (§ and obtain 
the density profile 


C{X) = iErfc(X) + E e-^V2Erfc(^-A^ (26) 


The density profile (261 holds everywhere and it does 
not depend on the initial density p. The latter prop¬ 
erty is known in the homogeneous case, so it must hold 
when X —>■ — 00 ; remarkably, it remains true even for the 
initially empty half-line. The dependence on the initial 
density can be observed e.g. for x ^ t, but for any fixed 
X this dependence disappears in the long time limit. 

On the interface separating the two half-lines, we have 
C{x = 0) = i -b ^ = 0.853553... (see Fig. [^. If there 

were no reactions, only diffusion, c{x,t) = ^Erfc^;^^, 


so on the interface separating the two half spaces the 
density is exactly a half that of the bulk. 

The exact solution ( |25[ ) contains information beyond 
the density. For instance, specializing ( [^ to E = 00 
gives the probability E{x;t) = E{x,y = 00 ;t) that the 
half-line {x, 00 ) is empty. This probability is a function 
of a single scaled variable: 



FIG. 3: The scaled density C{X) in all three processes and 
the scaled density of monomers Ci{X) for the aggregation 
process. The density C{X) is the same, Eq. (261, both for the 
coalescence and annihilation processes if the scaling forms are 
chosen according to ([^ for the coalescence process and (311 
for the annihilation process. In the aggregation process, the 
scaled densities of clusters of small mass are given by the same 
Eq. (|40[) as the density of monomers. 


The integral can be expressed through the error function: 

£(X) = l-lErfc(^^) (27) 

The behavior £(0) = ^ simply reflects that the right¬ 
most particle in the coalescence process is equally likely 
to be to the left and to the right of the origin [see also 

(pi)]. 


B. Irreversible Annihilation Process 


For the irreversible annihilation process ([^ , instead of 
empty interval probabilities one considers G{x,y\t), the 
probability that there is an even number of particles in 
the interval [x,y\. The probabilities G{x,y,t) satisfy the 
same equations (17)- as the empty interval proba¬ 
bilities (see m l35]). the only difference is the initial 


condition, viz. instead of (20) one gets 


Go(a^,J/) = ^ ^ X <( exp[2px] 


1 0 < X < y 

X < 0 < y 
exp[2p(a; — y)] x < y <0 


There is therefore a duality between annihilation and co¬ 
alescence, namely 


G{x,y,t) = \ + \E{x,y,t) 


(28) 


This duality holds if the initial densities differ by factor 
2 exactly: 


annihilation _ 1 coalescence 

P ~ 2 P 


(29) 


e{x-,t) = e{x) = 1 - 


2^/Tr 


due Erfc(M-l-A) 


This remarkable duality between annihilation and coales¬ 
cence extends to more complicated correlation functions. 
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It was found [3SH3S] in the context of homogeneous set¬ 
ting, and it also holds in our spatially-inhomogeneous 
setting. 

In the long-time limit, the initial densities are irrele¬ 
vant, to leading order, and therefore, we have 


= l-lErfc(i±J:)Erf(r^) 




ru^-Y 


du e 


dv e 




(30) 


The density has the scaling form 


1 




(31) 


with the same scaled density C{X) as before, Eq. (261. 

Specializing (30l to E = oo gives the probability 
£a(a;;t) = G(x,y = oo;t) that the half-line {x,oo) con¬ 
tains an even number of particles: 


£a(^) = l- 


4v^7_, 


due ^ Erfc(M-l-X) (32) 


In particular, the initially-empty half line contains an 
even number of particles with probability 

£a(0) = ^ Pa(21V) = ^ (33) 

N>0 


C. Irreversible Aggregation Process 


The average mass density = X)m>i 

is not affected by the aggregation process and hence it 
satisfies the diffusion equation dty = Dd^x^- Therefore 



To probe other average characteristics, we need more 
advanced methods. It turns out [331 131 ] that instead of 
empty interval probabilities it suffices to use Pm(x,y]t), 
the probability that the total mass inside the interval 
[a;, y] equals m. Assuming that all clusters are point par¬ 
ticles that diffuse with the same diffusion constant P, 
one finds that the probabilities Pm{x,y;t) satisfy 


The initial condition reads 

{ Sm,o 0 < X <y 

eP- x<0<y (38) 

My^^^pG-y) x<y <0 

The governing equations ( |3^ are linear, yet their anal¬ 
ysis is rather cumbersome, see Appendix Here we just 
present the major findings. In the scaling limit ([ 5 ]), the 
mass distribution acquires the scaling form 

Cm(Al, t) = Cm (AT) (39) 

For m = 0(1), the scaled densities become mass inde¬ 
pendent, i.e., the same as the monomer density 


C^(X) = Erfc (-^) + , 




(40) 


As a function of the scaled distance, the monomer density 
has a single peak (Fig.|^ and it is maximal at A « 0.84. 

One anticipates that the mass density Cm (A) depends 
on the scaled mass: 


0™(A) = $(M,A), M = 


(41) 


p\/27rPt 

We now use the sum rule X)m>iCm(A, t) = c{X,t), 
the scaled mass distribution (39), (41) together with the 
cluster density (§ to obtain an integral relation for the 
scaled mass density 

r °° p -^^/2 / s 

J dM$(M,A)=Erfc(A)-H ^ Erfc(^-^J (42) 

Similarly the sum rule Y,m>i''^'^rniX,t) = yi{X,t) to¬ 
gether with the exact result (34) for the mass density 
lead to another integral relation 

f fiMMd>(M,A) = Erfc(^) (43) 


The small mass tail of $(M, A) is given by (40); it would 
be interesting to derive the large mass tail. 


III. PARTICLE NUMBER 


d / 52 \ 

- 1 ) = C (^^ + ^ j »;<) 


The boundary condition ( |l^ generalizes to 
Pm{^: t) — 

The mass distribution is found from 

dPmix,y,t) 


Cmix,t) = 


dy 


y=x 


(35) 


(36) 


(37) 


Next, we study the total number of particles N{t) that 
infiltrate the initially empty half-line. This natural quan¬ 
tity is finite^ and it fluctuates throughout the evolution. 
Hence, the average does not fully characterize the proba¬ 
bility distribution function. The finiteness of N (t) follows 
from a heuristic argument. The decay law c ~ 
in the homogeneous case and the fact that particles infil¬ 
trate by distance of the order of V Dt into initially empty 
half-line suggest that N is indeed of the order of one. We 
can compute the average total number of particles (A) in 
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the initially empty half-line. For the coalescence process, 
we use Q to obtain 


poo -I poo 

(7V)c= / dxc{x,t) = — dXC{X) 
Jo Jo 


(44) 


which in conjunction with the density profile (26) leads 
to the announced result 0- For the annihilation process, 
the density is exactly two times smaller; accordingly, the 
average total number of particles is 


(iV)a = 4 + = 0.267077471... 

lb 47r 


(45) 


To learn more about the random quantity iV, one needs 
to compute its probability distribution P{N). The distri¬ 
bution functions for the reaction processes ([^ and ([^ are 
actually different. We are interested in both Pc{d^ and 
Pa{N). (We distinguish the two cases with subscripts.) 
The coalescence process is simpler, and a few exact re¬ 
sults are possible. 

The simplest exact result is the probability of finding 
no particles in the initially empty half-line: 


^c(O) = 


(46) 


To derive Eq. (46) consider the right-most particle R and 


notice that from its ‘view-point’ the coalescence process 
0 reduces to A + R ^ R. Thus the right-most particle 
is not affected by the coalescence process, so it performs 
a one-dimensional Brownian motion. With probability 
I the right-most is at a; < 0, and this is equivalent to 
saying that N = 0. 

The normalization requirement, ~ 


together with (46) indicates that J2n>i ~ I- 

Combining this result with {N)c = and 

Eq. 0 we obtain the sum rule 


y (iV - l)Pc{N) = 4 i = 0.0341549431... (47) 

• ^ /TT X 


N>2 


2tt 8 


This sum rule leads to the upper bound for the probabil¬ 
ity to find two particles: 


Pc(2) < 0.0341549431... 


(48) 


Combining the normalization requirement with (46) 
and (0 we similarly derive the sum rule 


Y,(N-2)PAN) = 1 + 1 

N>3 


16 47r 

arctan^;^^ — 2arctan("^ 


1 

= 0.00025091951... 


VVly. 


which yields the upper bound for the probability to find 
three particles in the initially empty half-line: 

Pc(3) < 0.00025091951... (49) 


A. Computation of Fc)!) 


We now consider the coalescence process and show that 
the probability Pc(l) to have exactly one particle in the 
initially empty half-line is given by (0. In principle, our 
procedure can be extended to Pc{N) with arbitrary fV; 
yet the computation quickly becomes prohibitive. 

To determine Fc(l) we use probabilities for two inter¬ 
vals to be empty. In Sec. |II A| we computed the prob¬ 
ability E(x,y) that the interval {x,y) is empty at time 
t. (In the following we do not display the time variable, 
so E{x,y) denotes E{x,y\t), etc.) Let E(xi,yi;x2,y2) 
be the probability that the intervals {xi,yi) and {x 2 ,y 2 ) 
are empty. We shall assume that xi < yi < X 2 < y 2 , 
so the intervals are non-overlapping. The probability 
E{1;2) = E{xi,yi;x 2 ,y 2 ) satisfies 

The solution of this equation, with obvious boundary 
conditions 

lim E{xi,yi-,X 2 ,y 2 ) = E{x 2 ,y 2 ) 
xltVl 

lim E{xi,yi;x 2 ,y 2 ) = E{xi,y 2 ) 

Vl1'X2 

lim E{xi,yi-,X 2 ,y 2 ) = E{xi,yi) 

X2ty2 

can be expressed through the single-interval empty prob¬ 
abilities [55] : 

E{1;2) = E{xi,yi)E{x2,y2) - E{xi,X2)E{yi,y2) 

+ E{xi,y2)E{yi,X2) (50) 


For our purposes, it suffices to consider a simple subset 
of empty interval probabilities, namely those with y 2 = 
oo. We denote E{x, y, z) = E(x, y, z, oo), in analogy with 


notation £(z) = E{z,oo) which we used in Sec. IIA 
Specifying (50) to this setting we get 


£(a:, y, z) = E{x, y)£{z)-E(x, z)£{y)+E{x)E{y, z) (51) 


Applying to Eq. ( j^ and taking the z ^ y limit, 
we find the probability R{x, y) that the two right-most 
particles are at x and y. 


R{x,y) 


dE{x,y) d£.{y) 


dx dy 


c{y) 


dlL{x) 

dx 


£( 2 /) 


d‘^E{x, y) 
dxdy 


(52) 


Further, the probability Pc(l) to have exactly one par¬ 
ticle in the initially empty half-line can be obtained by 
integrating the probability density R(x,y)\ 


poo pO 

Pc{l) = dy dxR{x,y) 
Jo J — oo 


(53) 


Computation of the integrals in Eq. (531, detailed in 
Appendix 10 leads to the the announced expression 0. 
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As a by-product of these calculations, we can deter¬ 
mine the probability density R{x, y) that the first and 
the second right-most particles are located at x and y. 
Using (52) we can express R{x,y) in the scaling form 


Rix,y\t) = 


1 




Ji{X, Y) = Erfc '^{X, Y) 


— e 


-X^l2 


Erfc(y) 


(54a) 


(54b) 


where we have used the shorthand notation 

4'(A:,y)= -F (y - a:)V^ ^ ^ 


l-iErfc(^) 


For the homogeneous coalescence process, the distri¬ 
bution of distance between adjacent particles is easy to 
determine, see e.g. nulls]. For the particles on the edge, 
however, the comp utati on is involved and the scaled dis¬ 
tance distribution (54b) apparently has not been known. 


For the annihilation process, even in the bulk (equiva¬ 
lently for the homogeneous setting) the distance distribu¬ 
tion has not been established despite of the considerable 
effort 


B. Simulation Results 

We performed numerical simulations to measure vari¬ 
ous statistical properties including in particular, those of 
the total number of particle N. In the simulations, we 
considered the lattice version of the diffusion-controlled 
coalesce and annihilation processes. Namely, we assumed 
that particles undergo (continuous time) random walk 
on the one-dimensional lattice, and instantaneously co¬ 
alesce (or annihilate) whenever two particles occupy the 
same site. In the long-time limit, the continuous (with 
particles undergoing Brownian motion) and the lattice 
versions lead to the same results. 

The initial condition become irrelevant in the long time 
limit, so we considered the simplest initial condition when 
half of the lattice is fully occupied. Thus initially every 
lattice site in the half-line a: < 0 is occupied and every 
lattice site in the half-line a: > 0 is empty. 


A 

Ua(A) 

Uc(A) 

0 

0.74 

0.50 

1 

0.25 

0.46 

2 

0.008 

0.03 

3 

1 

O 

CO 

3- 10-* 


TABLE I: The probabilities Pa(A') and Pc{N) for A < 3 as 
obtained from numerical simulations. 


We measured the stationary distributions Pa{N) and 
Pc(Af) for the number of particles in the initially-empty 
half space (x > 0) in the annihilation and coalescence 
processes. The results are listed in Table I. The numer¬ 
ical findings are consistent with the sum rule (33), viz. 



FIG. 4: Semi-logarithmic plots of the distributions Pa{N) and 
Pc(A) versus A^. 


^j,>gPa(2fc) = |. For the coalescence process, simula¬ 
tion results for Pc(0) and Pc(l) are in excellent agreement 
with theoretical predictions ( [46| and ([^. The simulation 
result for Pc(2) is close to the theoretical upper bound 
(48) which is natural since Pc(3) is very small according 
to another theoretical upper bound (49). Both distribu¬ 
tions Pa(A^) and Pc{N) apparently have Gaussian tails. 
The numerical evidence (see Fig. looks fairly convinc¬ 
ing, even though the range in N is very small. 


C. Characteristics of the leader 


Here we outline what we know about the position of 
the “leader” defined as the right-most particle. Features 
of interest include its average location {£), the probability 
density n(f, t) that the leader is at position £ at time t, 
etc. For the coalescence process 


(£)c = o, nc(f,t) 


1 

y/AnDt 


exp 


e - 

ADt 


(55) 


For the annihilation process, in contrast, {tjs. and t) 
are not known. Since the initially empty half-line is 
empty with probability exceeding ^, one anticipates that 
the average position of the leading particle recedes dif¬ 
fusively into the initially occupied half-line. This was 
indeed observed in simulations |32j where the probabil¬ 
ity density was also measured, and found to be 

asymmetric and clearly non-Gaussian |32j . 

We studied numerically the one-dimensional aggrega¬ 
tion process (§ with mass-independent diffusion coeffi¬ 
cients. We observed that the average mass of the leader 
(the right-most aggregate) grows as (m) = Bp^/Dt with 
P = 1.6 ± 0.1. This growth is anticipated, e.g., using 
(34) one finds that the average total mass in the initially 
empty half-line grows as p\/Dt/{A tt). 

The one-dimensional aggregation process generalizes 
the coalescence process, so it is more tractable than the 























annihilation process. The properties of the leader might 
be amenable to exact analysis. We already know the 
probability density for its location. It will be 

interesting to determine the mass distribution P{m,t) 
and the more detailed joint distribution i.e., 

the probability density that the leading particle has mass 
m and located at £. 


IV. DISCUSSION 

In summary, we studied three basic one-dimensional 
diffusion-controlled reaction processes—annihilation, co¬ 
alescence, and aggregation. In the case of aggrega¬ 
tion, we assumed that the diffusion coefficients are mass- 
independent, so it reduces to the coalescence process 
if one focuses only on the total cluster distribution. 
Our theoretical analysis generalizes the empty interval 
method to inhomogeneous initial conditions. 

We examined the evolution starting with the inhomo¬ 
geneous initial conhguration when a half-line is uniformly 
filled by particles, while the complementary half-line is 
empty. We computed the average particle density as a 
function of position. Our main finding is that while the 
overall density profile is a time-dependent quantity, sta¬ 
tistical properties of the total number of particles residing 
in the initially-empty half-space become independent of 
time, in the long-time limit. In particular, the total num¬ 
ber of particles in the initially empty half-line is finite and 
governed by a stationary probability distribution. We 
were able to compute analytically several features of this 
distribution function, e.g. the average. 

There are numerous basic questions about the behav¬ 
ior of one-dimensional diffusion-controlled annihilation, 
coalescence, and aggregation that remain unanswered. 
We limited ourselves to equal-time correlation functions. 
Many of them, e.g. {N{t)), become time-independent in 
the long time limit. Two-time correlation functions, such 
as {N{ti)N{t 2 )), are natural extensions of our study. 

Another interesting time-dependent quantities are var¬ 


ious first-passage (or persistence) characteristics [IHl IH] . 
For instance, the survival probabilities of the particles in 
the annihilation process decay algebraically with time, 
and the decay exponent depends only on the parity of 
the label corresponding to the initial order; these two 
exponents [32j haven’t been determined analytically. 

One can also ask about the probability S{t) that not 
a single particle ever entered the initially empty half¬ 
line during the time interval (0,t). For the coalescence 
process Sc{t) ^ since the problem reduces to the 

survival probability for a Brownian particle in one dimen¬ 
sion snun]. For the annihilation process Sn{t) ~ t 
as follows from equivalence with spin persistence problem 
[42] . The probability of never entering the initially empty 
half-line is equivalent to the probability that the position 
of the leader satisfies i{t') < 0 for all t' < t. Since the 
typical position of the leader scales diffusively, it is nat¬ 
ural to ask about the probability that i{f) < A-\/ Dt' for 
all t' < t, where A is some fixed constant. This probabil¬ 
ity decays algebraically, S{t,A) ^ as t ^ oo. For 

the coalescence process we have effectively a single parti¬ 
cle problem which is solvable and the exponent 9c{A) is 
known |43j for arbitrary A; for the annihilation process, 
only 0a(O) = ^ is known [35]. Another natural gen¬ 
eralization involves first-passage properties involving the 
number of particles residing in the initially-empty half 
space [44] . 

The lack of analytical tools for probing the properties 
of the annihilation process is frustrating, although not 
surprising in the light of earlier work [571 - I55| . We think 
that establishing analytical tools for probing the anni¬ 
hilation process represent the most important and most 
challenging extension of the current work. Some deli¬ 
cate properties of the homogenous annihilation process 
have been derived analytically |42| using field-theoretical 
methods; perhaps, the techniques of [42] can be extended 
to the inhomogeneous setting. 

We acknowledge support from US-DOE grant DE-AC52- 
06NA25396 (EB). 
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Appendix A: Scaling Behavior in the 
One-Dimensional Aggregation Process 


One can try to solve (35) using the generating func¬ 
tion technique. Here we use another approach relying on 
the recurrent nature of Eqs. (35). We first show in detail 
how to find the density of monomers and then generalize. 
The probability Pi(x,y,t) that the total mass contained 
in the interval [x, y\ is equal to one is clearly the proba¬ 
bility that the interval [a;, y\ contains one monomer. This 


probability varies according to diffusion equation 


dt 


P,=D 


dx^ 


dy'^ 


P^ 


supplemented by the boundary condition 
Pi{x,x;t) = 0 


(Al) 


(A2) 


As earlier, we extend the definition of Piix, y, t) from the 
physically relevant half-plane x < y to the entire plane. 
The initial condition 


Pi(t = 0) = 


0 

{—px) 

py 

p{y — x) 
—p{x — y) 


X > 0, y > 0 
X < 0 < y 
y < 0 < X 
X < y < 0 
y < X < 0 


(A3) 


agrees with (38) when x < y, while in the supplementary 

assures that 


half-plane x > y the choice made in Eq. 
the boundary condition (A21 is manifestly obeyed. 


Solving the diffu sion equation ( |Al[ ) subject to the ini¬ 
tial condition (A3) yields 

A-nDt 7 ° 

-Pi = dxj dyo{-xo)eP^«G 

P J— oo Jo 

pOO pO 

-F dxo dyoyoe^^°G 

Jo J — oo 

/ O pyo 

dyo / dxo {yo - xq) 

-OO J —OO 

/ O pXQ 

dxo / dyo {yo - xo) 

-OO J —oo 

Here G = G{x,y,t\xo,yo) is the diffusion propagator: 

{x - xo)^ + (y - yo)^ 


G{x,y,t\xo,yo) = exp<^ - 


ADt 


Recalling (37) and performing straightforward calcula- 
)non 

dyo 


tions one finds the monomer density ci{x,t) 

8TT{Dt)'^ 


-Cl = 


dww'^e-P^Hi (A4) 

Jo Jo 

pOO pOO 

/ dyo dxoXoe~P^°{yo- xo)H 2 
Jo Jo 


where we have used shorthand notation: 

Hi = H{x\- yo,-yo - w) 

H2 = H{x\ - xo,yo) 

(x - Xq )^ + (x- yoY 
ADt 


H{x\xo,yo) = exp<^ - 


The integral representation (A4) is exact, that is, valid 
at all t > 0. In the long-time limit, the monomer density 
simplifies. Let us first compute ci(0,t), the monomer 
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density exactly on the interface. When a; = 0, the right- 
hand side of (A4) becomes 


Vo + iVo + 


POO POO r 

J dyoj exp| — 

POO POO r 

+ J dyo J dxQXQe~P'^°{yQ- xo)eyi^^- 


^l + vl 

ADt 


In the long-time limit, more precisely when p^Dt 3> 1, 
the integral in the second line dominates and asymptot¬ 
ically it grows as 2Dt/p^. Therefore 

1 


ci(0,f) = 


AirpDt 

More generally, ci{x,t) = {A'KpDt)~^/'^ for x ^ p~^. 

In the scaling limit ([^ with X being fi xed, the integral 
in the first line on the right-hand side of (A4) approaches 
to Erfc(X), while the integral in the second 

line of (A4) tends to 


2Dt 

„2 






Thus 

Ci{X,t) = 


1 

AnpDt 

1 




1 


, „ , ,_=Erfc(X) (A6) 

The term in the second line is negligible for finite X, 
both positive and negative, but we kept this term as it 
prevails when X <C — \J\a[p‘^Dt). This term determines 
the asymptotic monomer density far away from the inter¬ 
face (equivalently, for the uniform initial distribution): 

ci(—oo, t) = p~^ 

A similar treatment leads to 


(A7) 


8Tr{Dty 

pm 


dyo 


dw 


,,m+l 


ml 


“iji 


-k 


/ dxo / 
/o Jo 


dyo -X e 
ml 


0 _ ^-pXQ 


{Vo - Xo)H2 


For clusters of finite mass, m = 0(1), one gets 

Cm ( 0 — 


AirpDt 
m + 1 




Ixe-V/^Erfc(-X)+e 

The term in the first line dominates for finite X, and it 
does not depend on the mass in the leading order. 


Appendix B: Computation of the integrals in (531 


Plugging (52) into (531 and performing the integration 


over X we get 

POO 

Pc{l) = / 

^0 

poo Q 

- J dy8{y) — [£'( 0 , y) - E{-oo, y)] 


dy [£(0, y) - £(-oo, y)] 


Since E{—oo,y) = £(—oo) = 0 and £(0) = h we have 


pOO 

Pc(l) = dy 

Jo 




dy 
dy c{y) 


dy 


(B2) 


Integrating by parts, we simplify (B2) to 




Specifying (251 to a: = 0, we get 


£(0,2/) = l-iErfc(f)Erf(|) 

^ poo ^ pu-\-Y 

TT }_X. 


/ oo pli 


Expressing the integral through error functions yields 


£(0,2/) = 1 - iErf (f) + i[Erf (f )]^ - ^Erf (^) 


It is convenient to transform y to Y in the first integral 
in (B3). Equation (27) gives dy = dY 




The above results togethe r wi th Jp dyc{y) = (£)c, see 
(44), allow us to re-write (B3) as 


^c(i) = 7- 


^OO -I 

dze-2^'n(z)- -(iV)c 


(B4) 


with 


n(z) = 1 - iErf(z) + 7Erf(z)]" - \Eri{V2 z) 


Computing the integrals in (B4) one gets 


dze =1 


*e-“’Erf(..) = Mtct.„(Aj 

[ dz Erf(y2 z) = ^ 

V TT Jo 2 

'Hi |Erf(*)l^ = ?™ta„(i) 


Using these expressions and recalling 0 we reduce ( |B4[ ) 
to the announced result d^. 


dyc{y)[8{0) - £(-oo)] 


(Bl) 










































